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O ' One of the greatest challenges in cosmology today is to determine the nature of dark 

^ , energy (DE), the source of the observed present acceleration of the Universe. Besides 

■ the vacuum energy, various DE models have been suggested. The tests that have been 

proposed to differentiate among these models are based on observations of galaxies at 
high redshift (z > 0), to be obtained in the future. We suggest here a new test that is 
valid at z ~ 0. It is based on existing observational data, numerical simulations, and three 
well known analytic models that evaluate the bias parameter 6, the ratio of galaxy to 
' dark matter (DM) fluctuations. These analytic models are based on the physical processes 

' involved in the formation of stars and in the formation and merging of galaxies. The value 

of b(z) obtained in each model is a function of the DM growth factor D(z), which, in turn, 
is a function of the DE. We show that the equations for b in all three analytic models can 
be reduced to the form of a known constant plus the term E[D(z = 0) / D(z)] a , where 
a = 1 or 2 and E is a free parameter. Using the value of b obtained by the 2dFGRS 
consortium for the ACDM model, to normalize E, we find that all three analytic models 
predict 6 2 (0) = 1±0.1 for all DE models. Since we use the result that b 2 (0) ~ 1 from the 
2dFGRS consortium for the ACDM model, the L» galaxy used in our text is the same as 
that obtained by the consortium for a broad range of galaxy types, using the Schechter 
function fit to the overall luminosity function of their ~ 220, 000 galaxies. Numerical 
simulations that evaluated b 2 (0) for the ACDM and CDM (A = 0) models also obtained 
6 2 (0) = 1 ± 0.1. Since this value of b 2 (0) is indicated by numerical simulations as well as 
by all three popular analytic models, which are normalized by the 2dFGRS consortium 
result for the ACDM model, we suggest the condition that b 2 (0) = 1 ± 0.1 at z = as 
a new test for the viability of DE models. Thus, for a given observed galaxy fluctuation 
spectrum such as that of the 2dFGR,S consortium, if the DM fluctuations are greater or 
less than the galaxy fluctuations by more than 10%, the DE model can be discarded. 
As examples of this test, we show that three popular DE models do not satisfy this 
test: the vacuum metamorphosis model deviates from b 2 (0) = 1.0 at z = by 20%, the 
brane-world model by 26% and the supergravity (SUGRA) model by 38%. 

Keywords: Observational cosmology; Dark energy; Spatial distribution of galaxies. 
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1. Introduction 

The nature of dark energy (DE), the existence of which was first indicated to explain 
the recent SNIa observ atio n of the acceleration of the Universe, is one of the major 
problems in cosmology. Theories in which gravity is modified as well as those that 
include parametrizations of the DE equation of state (EOS), w(z) = p/p, where pip) 
is the pressure (energy density) of the DE, have been suggested to explain it.^HSI 
Based on observations, various constraints have been put on the EOS for a variety 
of models (see e.g., Refs. I7HT3|). In order to investigate DE models, we used the bias 
parameter b 2 , the ratio of galaxy clustering £ gg to dark matter (DM) clustering £,dm- 

Observed galaxy clustering £ gff [or its power spectrum P(k)] does not directly 
provide information about the DE. The DE affects the DM clustering £dm which, 
in turn, affect ^ gg . In order to obtain information about DE, we need to know 
the ratio b 2 = £, gg /£,dm- In the future, observations will be used to determine £,dm- 
Thus, along with observations of £ 9S , we will be able to obtain b 2 . One very impor- 
tant observational project to determine ^ m is the DES (the Dark Energy Survey) 
(http://www.darkenergysurvey.org). DES will study the growth of DM fluctuations 
as a function of redshift with "weak" gravitational lensing, produced by the DM 
fluctuations. The results of DES will take some time, however, to become available. 
In the present paper, we use observations to normalize existing analytic models to 
predict b 2 (0) ~ 1 at z = 0, which we propose as a new test for DE models. At 
the present time, this test is very useful for limiting viable DE models. Recently, 
Grande et al. used it to limit the validity of two DE models. ^31 

Somerville et al. define the general bias parameter as b(6)6 =< S g \5 >= 
J d5 g P{8 g \5) 8 g for a DM fluctuation <5. 021 The bias parameter b{8) is the average of 
the probability P(8 g \S) that there is a galaxy fluctuation S g within the matter fluc- 
tuation 5. This relation fully characterizes the mean non-linear biasing and reduces 
to the linear biasing relation, b5(X) = 8 g (X), if b is independent of S, where S g (X) 
is the galaxy fluctuation in a sphere of radius A. In order to track the formation 
of galaxies and quasars in their simulations to evaluate 6, Somerville et al. used 
a semi-analytic model to follow gas, star, and supermassive black hole processes 
within the merger trees of DM halos and substructures. This semi-analytic model 
is described in Refs. [TBI H3 HH1 The modelling assumptions and parameters were 
adjusted in order to fit the observed properties of low redshift galaxies, primarily 
their joint luminosity-color distribution and their distributions of morphology, gas 
content, and central black hole mass. 

The observed galaxy power spectrum from the final 2dFGRS catalogue can be 
found in Ref.Qj)] (Cole et al.). According to Cole et al., the large-scale linear bias 
factor for L* galaxies, k < 0.1/iMpc -1 is b = 1.03 for the ACDM model. This 
is consistent with their previous result for brighter L$ galaxies, b(L$, z = 0) = 
1.10 ± 0.08, obtained from APM-selected massive galaxies (L s = 1.9L*).^ 

Bias depends on the luminosities of the galaxies studied and is known as "lu- 
minosity segregation" . In order to normalize the bias to L* galaxies, the 2dFGRS 
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consortium used the expression b(L) = 0.85 + 0.15(L/L») to take into account 
luminosity segregation.^^! Since L* galaxies are almost unbiased, the 2dFGRS con- 
sortium result indicates that b 2 (z = 0) ~ 1.0 to a 10% accuracy for the ACDM 
model. We use this important result in our new test for DE models. 

The comparison between the DM fluctuations at the CMB epoch and those at 
the present, involves two steps: 

(i) the calculation of the DM perturbation amplitude at the CMB last scattering 
epoch ^dmis from the observed Ci of the CMB and the cosmological parameters; 
and 

(ii) the calculation of the DM perturbation amplitude at the present epoch ^dmo 
from the growth rate and £,dmis- 

We assume negligible DE at the recombination era in the present paper, making 
the perturbation amplitude at recombination independent of DE models. This is a 
reasonable assumption since the flatness of the Universe, which we assume for all 
our DE models, primarily determines the sound horizon for the CMB data and the 
mapping between the matter power spectrum k and the CMB I's [step (i)], while 
the DE primarily determines the growth rate [step (ii)] . 

The perturbation amplitude at the recombination era depends primarily on the 
cosmological parameters f^Mi ^b> and H in a flat Universe. The parameter 
is generally attributed to neutralinos produced in the supersymmetric dominated 
primordial Universe, is produced in the baryogenesis primordial Universe, and 
the present Hubble parameter Hq in a flat Universe is due to £l M and £l B . 

These parameters can be obtained from observations, using methods that are 
not in any way at all connected with the CMB and DE models. The Hm can 
be obtained from high precision observations of galaxy clusters, using the relation 
= (M/L)J/p c> where M(L) is the mass (luminosity) of a given cluster, J is 
the luminosity density in the field around the cluster, and p c is the critical density. 
Other methods using clusters and galaxies that depend on the linear theory growth 
factor (and thus, on the DE model), are also used to determine f2j\/. The parameter 
VLm can be obtained from the observed number of clusters as a function of redshift, 
compared with the number predicted by the Press-Schechter relation or the Shcth- 
Tormen model, which depends on the linear theory growth factor and the DE 
model. Yet another method of obtaining Qm, which depends on the linear theory 
growth factor, uses as, the present average amplitude of the dark matter fluctuations 
in a sphere of radius 8/i _1 Mpc. The parameter tts can be determined from high 
precision deuterium abundance observations of quasar absorption lines, and the 
present parameter Hq, from data of Cepheids in nearby galaxies. 

Analytic models that evaluate the linear bias evolution b(z) as a function of 
redshift, depend on the EOS of the DE.^21 These analytic models predict b 2 (z ~ 
0) ~ 1.0 for the flat DE ACDM model. The reason for this was noted by Tegmark 
and Peebles:^2l "...even if galaxies initially were uncorrelated with the mass, they 
would gradually become correlated as gravity draws them toward overdense regions 
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and one might expect this process to drive b toward unity". Numerical simulations, 
such as those of Somerville et al., also obtain this value for 6 2 (0) for very different 
DE models, for example, the CDM (A = 0) and ACDM models. 

From the analysis of the 2dFGRS data by Cole et al., b\ CDM (0) ~ 1.06. Nor- 
malizing the analytic bias models, discussed in Section 2, to 0.9 < b\ CDM (0) < 1.1, 
we find that 0.9 < 6 2 (0) < 1.1 for viable DE models, which we propose as a new test 
for DE models. Since we use the result that b 2 (0) ~ 1 from the 2dFGRS consortium 
for the ACDM model, the L* galaxy used in our text is the same as that obtained 
by the consortium for a broad range of galaxy types, using the Schechter function 
fit to the overall luminosity function of their ~ 220, 000 galaxies. 

We apply our DE model test to several popular DE models. We relate the b 2 of 
the DE model at z ~ to a factor, F = \b 2 — b 2 A )/b 2 A \, where b\ is the b 2 of the 
flat ACDM model at z ~ 0. A 10% deviation of b 2 from b\ normalized to unity 
implies a 10% deviation of F from zero, F = 0.1. Viable DE models then need to 
have F < 0.1 (i.e., b 2 = 1.0 to a 10% accuracy at z = 0). 

Three well-known analytic models of linear bias evolution are discussed in Sec- 
tion 2. We discuss the effect of DE on the linear growth of (Sp/p) in Section 3. 
DE models can be described by an EOS, P/p = w(a) = wq + w a (l — a), where 
P{p) is the pressure (energy density) at the cosmic scale factor a)** It is shown 
in this section, that the permissable values of F < 0.1 limits the parameters uuq 
and w a . We then go on to disc uss three popular DE models: the fivc-dimcnsional 
brane- world model (BWM),^31the vacuum metamorphosis model ( VMM), ^1 and 
the supergravity (SUGRA) model. Conclusions and discussion are presented in 
Section 4. 

2. Analytic bias models as a function of DE 

As noted above, the 2dFGRS final results as well as the numerical calculations of 
Cole et al. indicate a bias parameter b{z = 0) ~ 1.03 for the ACDM model. From 
these results, we have b\ CDM ~ 1.0 to better than a 10% accuracy. We use this 
result to normalize well-known analytic models for b 2 {z) in order to obtain 6 2 (0) for 
viable DE models. 

2.1. Galaxy Conserving Model (GCM) 

In this model, galaxies behave as test particles, with their intrinsic properties con- 
served (see e.g., Refs. [23J [27l [28)l . An L* galaxy is a present massive galaxy. It 
probably started to form as a small galaxy at a high redshift, gradually building up 
by the process of merging. The assumption of the creation of an L* galaxy at a high 
redshift, with no further creation or merging, is an approximation of the evolution 
of galaxies. The linear bias parameter relating the density distribution to its 
mass density for a galaxy population formed at a given cosmic epoch is given by 



bacM(z) = 1 + [&(**) - 1 



D(z) 



1 + (&o - 1) 



D(z = 0) 
D(z) 



(1) 
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where bo is the bias parameter at the present epoch and D(z) = Sp/p is the linear 
growth of the density fluctuations (details of D(z) are presented in the Section 3). 

2.2. Galaxy Merging Model (GMM) 

The evolution of galaxy clustering is associated with host dark matter halos in 
this model. We use the analytical expression obtained by Ref . [29] for the halo- halo 
correlation, 

Z hh (r,M) =b 2 (M)Um(r), (2) 

loni 

where the bias parameter, according to the Press-Schechter formalism, 1 ^ can be 
written as 

1 n / 2/r s D 2 (z = 0) 

where v is the ratio of S c to the average value of the DM fluctuation a(M), for 
the mass M, [Eq.(5)]. From the dynamics of the spherical collapse in an expanding 
background, the factor 6 C was derived, and shown to be S c ~ 1.69, independent of 
the DE model by Ref. EU The factor 

u 2 /S c = 5 c /a 2 (M) (4) 

is independent of z and a 2 (M) is defined as 

(M) = D 2 (z)/D 2 (z = 0) ' (5) 

where 

a i( M , z) = D2 wivy = °) r k 2 p(k)w 2 (k, M)dk (6) 

2?r Jo 

is the rms of the linear density fluctuation of top-hat spheres containing an aver- 
age mass M. The function P{k) is the linear power spectrum at redshift zero and 
W(k, M) is the Fourier space representation of the real space top- hat enclosing the 
average mass M. 

2.3. Star Forming Model (SFM) 

DM halos with masses greater than a given mass M can be identified with galaxies 
with luminosities greater than a corresponding luminosity I at a redshift z in the 
star-forming model (see e.g., Ref. 32). In this model, the bias evolution is 

D 2 (z = 0) 
D 2 (z) 

The equations for b(z) in the above three models (i.e., Eqs. (1), (3) and (7)) can 
all be written in the form 

~D(z = 0) 



SFM (M, z) = l+ y/S c ) . (7) 



b(z) = A + E 



D(z) 



(8) 
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where A is a known constant, a = 1 or 2, and E is a free parameter. For the ACDM 
model obtained by the 2dFGRS consortium, b 2 = 1.0 ± 0.1 at a mean redshift 
z = 0.17.HSI xhe curves for the star forming and galaxy conserving models in Fig. 1 
were normalized to 6 2 (0) = 1.1, the approximate maximum, value of the 2dFGRS 
consortium. From Eqs. (1), (7) and Fig. 1 we note that the galaxy conserving model 
has b 2 equal to a constant plus a term cx z and the star forming model a constant 
plus a term cx z 2 . Both models predict a change of b 2 from z — 0.17 to z — by a 
very small factor ~ 0.1%. The 2dFGRS result b 2 = 1.0 ± 0.1 can then be used at 
z = 0: 



6 2 (0) = 1.0 ±0.1. 



(9) 



We use Eq. (9) to normalize E in Eq. (8) for all DE models. 

This article is based on the requirement that b 2 (z — 0) = 1 to an accuracy of 10% 
for all viable DE models, which is consistent with the results of the above analytic 
models. It is also consistent with the numerical calculations of Somerville et al., who 
showed that b 2 ~ 1.0 for both the CDM and ACDM models, (fi M = 1-0, = 0) 
and (Qm = 0.3, f^A = 0.7), respectively. 

We can compare the predictions of the above analytic models against observa- 
tions of galaxy formation as a function of redshift. In Fig. 1, we show the galaxy 
conserving model (short-dashed curve) and the star forming model (continuous 
curve) for ACDM, normalized to b 2 = 1.1 at z = 0. The curves are in agreement 
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Fig. 1. Analytic models of linear bias evolution, b 2 (z) for the fiat ACDM model (CIm = 0.3, Qa = 
0.7), normalized at & A (z = 0) = 1 with a deviation F = +10%, compared with the b A for the 
volume-limited sample data with -Mg = —20 + 51og/i from the first-epoch VIMOS-VLT Deep 
Survey (VVDS). The solid and the short-dashed lines correspond to the star-forming and galaxy 
conserving models of b 2 (z), respectively. The long-dashed line is b 2 (z) = 1 for comparison. The bias 
at z = 0.17 (the effective depth of the 2dFGRS survey) was inferred from the 2dFGRS galaxies,E31 
which had the same median luminosity of the volume-limited VVDS sample, L/L* ~ 2. 
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with the ACDM b 2 (z) from the first-epoch VIMOS-VLT Deep Survey (VVDS) in 
Ref.[32l (Marinoni et al.). (We do not plot the merging model since, as was shown 
by Marinoni et al. , it does not describe the redshift evolution of the bias of their 
galaxies well.) 

In Fig. 1, the values of E and A in Eq. (8) are E = 0.1 and A = 1.0 for 
the galaxy conserving and star forming models. Using the value of E for the star- 
forming model in Eqs. (4)-(7), we find that the dark matter fluctuations for the 
massive galaxies, a(M) =4.1, its present linear extrapolated value. The fact 
that a(M) = 4.1 is greater than 6 C = 1.69, implies that the halo of the L* galaxy 
has already collapsed, which is reasonable. 



3. DE and the growth of density fluctuations 

The nature of DE is still unknown and there are many alternative models which try 
to explain it. In addition to the popular cosmological model of a constant vacuum 
energy, described by a cosmological constant, there are models that modify gravity 
as well as DE those that parametrize the DE EOS, w(a) — wo + w a (l — a), setting 
values for wq and w a . For these models, the Fri edmann equation can be written in 
a general form in terms of an effective EOS. 1 ^ Modelling the DE as an ideal fluid 
in a flat Universe, we can write the Friedmann equation as 

= a° M {i + z f + (i - n°,) e 3J *dWi+*')ii+M*')] } (io) 



H 2 



or 



Hi 



5H 2 
Hi 



(11) 



where Hq is the present value for the Hubble parameter, fl M is the present normal- 
ized matter density, and SH 2 /Hl depends upon the DE model. The EOS w(z) for 
the DE can be written as 



, ld\n5H 2 /H2 



(12) 



Linear growth of a density fluctuation, D = Sp/p, depends on the EOS. We define 
the growth factor, G = D/a, where the cosmological scale factor is a = 1/(1 + z) 
and G is normalized to unity at z ~ 1100, the recombination epoch. ^ In terms of 
G, we have 

G'(a) 3 1 - 10(0) G{a) 



G"(a) + 



3 w(a) 



2 1+X(a) 



2 1+X(a) 



0. 



(13) 



where X(a) is defined as 



X(a) 



M • 



SH 2 /H 2 



For large X we recover the matter dominated behavior D 



(14) 
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It is generally required that the DE was very small compared to the cold DM 
for redshifts z > 10 in order for the latter to create the structure in the Universe 
at redshifts z < 10. In this paper, we analyze only dark energy models in which the 
DE density at the recombination redshift or higher (z > 1100) was negligible, as in 
the ACDM model. 

One of the DE models studied was the Vacuum Metamorphosis Model (VMM), 
which has zero DE at recombination. A second model was the Brane World Model 
(BWM), in which the ratio of the DE density to the DM density is proportional 
to (1 + z)~ 3 and negligible at high redshifts. A third model was the supergravity 
SUGRA model, which becomes a CDM model at high redshifts with negligible DE 
at recombination. 

Much effort has been made to obtain the density power spectrum normalization 

at the recombination era. For example, Spergel et al. used the three year WMAP 

data with a ACDM model to obtain the density power spectrum normalization As 
1041 

at recombination. 1 ^ 

Our DE model test is independent of As since it is primarily dependent on the 
evolution of the DM and radiation densities at z = 1100, the recombination era, or 
higher, when DE was small or negligible. The models that we examine here have 
negligible DE for z > 1100 and, thus, a negligible effect on Ag. It is to be noted that 
the ACDM model, used by Spergel et al., is one such model, in which DE density 
was one billionth that of the total energy density at the recombination era. 

For our DE model test, it is not necessary to obtain the exact value of the 
normalization factor. At a given epoch, the only relevant difference between the 
different models studied is the value of the DE density, all other factors being equal. 
Since at the recombination epoch, all of the models we examined had a negligible 
DE density, they are essentially identical and their effect on the normalization factor 
is identically negligible. Thus, an uncertainty in the normalization factor (which is 
not known to better than ~ 20%) does not enter into the uncertainty of ~ 10% in 
the factor F [Eqs. (15) and (16)], used in our DE model test, since all models are 
normalized at the recombination era. 

It is important to emphasize the difference between the objective of Spergel et 
al. and ours. The objective of Spergel et al. was to obtain a value for As as accurate 
as possible as well as other parameters, such as the normalized DM density Qm, 
using one DE model, the ACDM model. Our objective here, is to discard DE models 
which have negligible DE densities at the recombination era, but that dominate the 
total energy density of the Universe at present. 

Galaxies began to form at a redshift z ~ 20, growing until the present time 
(z ~ 0). The formation and growth of the galaxies, described by their fluctuations 
£gg, are greater for larger DM fluctuations ^ m in the redshift interval z ~ 20 to 
z ~ 0. Since ^ m monotonically increases with a decrease in redshift and DE does 
not effect it at very high redshifts (z 3> 20), £d m (0) is sensitive to the DE model 
and the growth factor. Our test, b 2 (0) ~ 1, does not tell us whether the growth 
factor or (0) needs to be big or small. What b 2 (0) ~ 1 does say is that whatever 
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£ gg (0) is, £dm(0) must accompany it such that 6 2 (0) = £gg(0)/£dm(0) 
We define the deviation from the standard ACDM model by 



F 



D 2 



r> 2 



D 2 



z=0 



(15) 



where D\ are the density fluctuations in the standard flat ACDM model and D 2 
are the density fluctuations of an arbitrary DE model. The functions D/a and 
D\/a are normalized to unity at the recombination era, z ~ 1100, where a = 
1/(1 + z) is the cosmic scale factor. Therefore F in Eq. (15) is concerned only with 
the relative growth of the density fluctuations between a given DE model and the 
flat ACDM model, from the recombination epoch to the present era. It thus has the 
important characteristic of being independent of the power spectrum normalization 
at recombination. 

The linear bias parameter is defined by b 2 = P gg /P mm , where P gg is a galaxy dis- 
tribution power spectrum (e.g., that of 2DFGRS or SDSS) and P mm {z) oc (5p/p) 2 = 
D 2 {z) is the DM power spectrum. The difference in the DE models is in the evolu- 
tion of Sp/p, which occurs after the recombination era at redshifts z < 10. For this 
reason, we use the factor F, which examines the effect of DE models on the evolu- 
tion of Sp/p after the recombination era. The DE models that we investigate do not 
effect the evolution of Sp/p before the recombination era and, therefore, (Sp/p) rec 
is the same for all the DE models studied. Thus, Eq. (15) becomes 



F 



z=0 



(16) 



A maximum 10% deviation of b 2 from b\ with b\ 



1.0 at z = implies that the 



maximum F is F„ 



0.1. 



Table 1. Best fit values of wo and w a 
for the EOS, w(a) = i»n + w n (1 — a), 
for the Gold SNIa dataset ES w jth a de- 
viation F = 0.10 ± 0.02 and a matter 
density Q.° M = 0.28 ± 0.02. 



W a 


w 


1.53 


-1.72 ±0.02 


1.63 


-1.77 ±0.02 


1.73 


-1.82 ±0.02 


1.83 


-1.86 ±0.02 


1.93 


-1.82 - 0.01 


2.03 


-1.89 - 0.01 


2.9 


-1.86 - 0.02 
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3.1. DE models described by a parametrized EOS 

We first discuss a parametrization for the EOS, w(a) = wo + w a (l — a), which has 
been widely used for DE models since it is well-behaved at high redshifts, unlike 
w(z) = wq + wiz, for example, which diverges at high z. This parametrization was 
introduced by Ref. The best fit parameters, Wq and w ai which are consistent with 
the Gold SNIa dataset, were found to be in the intervals —1.91 ^ wq ^ —1.25 and 
1.53 ^to„< 5.05.ES Assuming that F = 0.10 ± 0.02 and D,° M = 0.28 ± 0.02, we 
obtain the best fit values, —1.91 < wq < —1.72 and 1.53 < w a < 2.9, shown in 
Table [TJ 

Table 2. The deviation F for the BWM and 
the VMM, respectively, as a function of the 
matter density fi)L. The values Hqt c , m? and 
Zj are defined in § 3.2. 





Fbwm 


Hor c 


-FVmm 


lu- 


z j 


0.26 


0.27 


1.4 


0.29 


ll 


1.4 


0.28 


0.27 


1.4 


0.24 


11 


1.4 


0.3 


0.26 


1.4 


0.20 


11 


1.3 


0.32 


0.25 


1.5 


0.16 


11 


1.2 


0.34 


0.25 


1.5 


0.13 


11 


1.2 


0.36 


0.24 


1.6 


0.095 


10 


1.1 


0.72 


0.11 


3.6 


0.24 


8 


0.5 



3.2. Brane-world and vacuum metamorphosis models 

In the BWM,I231 gravity is modified by adding a five-dimensional Einstein-Hilbert 
action, that dominates at distances that are larger than the crossover length r c , 
which defines an effective energy density, flbw = (1 — ^m) 2 /4 = l/(4^o r c)) ^ or a 
flat Universe. The factor SH 2 /Hq in Eq. (fTTj) then becomes 

5H 2 /Hl = 2£l bw + 2^Z\/n° M (l + z)3 + n bw . (17) 

In the VMM,I2S1 the vacuum contributions are due to a quantized massive scalar 
field, which is coupled to gravity. For z < Zj, 5H 2 /Hq in Eq. (fTTj) is 

SH 2 /H 2 = (1 - m 2 /12)(l + zf + m 2 /12 - tt° M (l + zf , (18) 

where z 3 = [m 2 /(3fl a M )] 1 / 3 - 1 and m 2 = 3fl Q M [(A/m 2 ) - (1/3)]- 3 / 4 . Both the 
BWM and the VMM can be described by the EOS, w(a) = Wp + w a (l - a), with 
(w Q ,w a ) = (—0.78,0.32) and (w ,w a ) = (—1,-3), respectively.^ 

It can be seen from Table 2 that an agreement within 10% between the VMM 
and the ACDM model is possible only for a matter density VL° M « 0.36. For the 
BWM, an agreement with the ACDM model within 10% is possible only if the 
matter density Vl a M w 0.72. However, both of these values for f2° f are greater than 
the observed value. 
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3.3. Supergravity model 

The SUGRA model is very attractive for explaining the acceleration of the Uni- 
verse.^It can be described by the EOS of § 3.1, with w a = -0.82 and w a = 0.58.E3 
This equation of state is in agreeme nt wit h observations for the low redshift SNIa 
dataset and galaxy distribution data.fHSHI However, F for this model makes it invi- 
able. 

Figure 2 shows the growth of the density fluctuations G — [(Sp/p) /a] as a 
function of a for the BWM, the VMM, and the ACDM model. It shows that the 
growth of Sp/p is smaller for the BWM and SUGRA model than for the ACDM 
model. The deviation F for the SUGRA model is F ST jgrA « 0.38 for Q.% = 
0.30Zq'q2, which gives a bias parameter, normalized to the b\ CDM (z = 0) = 1, 
appreciably greater than unity: b 2 SUGRA (z = 0) = 1.38. 

4. Conclusions and Discussion 

In this paper we suggested a new test for the viability of DE models, based on 
the value of the bias parameter b, the ratio of galaxy to DM fluctuations, at z = 
0. If it were the case that we knew nothing about galaxy formation, b could, in 
principle, have been anything at all, i.e., very much less or very much greater than 
unity. However, our present knowledge of galaxy formation from analytic models and 
numerical simulations, indicates that b 2 is close to unity at z ~ 0. This information 
can be used to discard DE models that do not predict b close to unity at z = 0. 

We studied three popular analytic models for b(z). We showed that the equations 
for b in all three analytic models can be reduced to the form of a known constant 
plus the term E[D(z = 0)/ D(z)] a , where a = 1 or 2 and E is a free parameter. 




0.2 0.4 0.6 0.8 1 



a 

Fig. 2. The growth of density fluctuations, G = [(Sp/p) /a], for the VMM, ACDM, BWM, and 
the SUGRA models. The dashed lines show the deviation of the matter density, = 0.28±0.02. 
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Using the value for b obtained by the 2dFGRS consortiums^ for all three models 
predict b 2 (Q) — 1.0 ± 0.1 for all DE models. This value of b is also in agreement 
with numerical simulations that evaluated b 2 (0) for the ACDM and CDM (A = 0) 
models.^ Since this value of b 2 (0) is indicated by numerical simulations as well 
as by all three popular analytic models, which are normalized by the 2dFGRS 
consortium result for the ACDM model, we suggest the condition that b 2 (0) — 1±0.1 
at z = as a new test for the viability of DE models. 

Obtaining 6(0) from galaxy observations involves a complex process of combining 
data from all types of galaxies (see e.g., Ref. 39). These complexities reflect the 
galaxy formation process. As in the standard analysis for the evaluation of the bias 
parameter b (see e.g., Cole et al. Ref. 19}, we assume that b is independent of scale 
for k ^ 0.1/iMpc -1 and that all galaxies are normalized to a standard massive 
bright galaxy with luminosity L*. Cole et al. made the following normalization of 
the bias parameter to a galaxy of luminosity L* for the galaxies in the 2dF Galaxy 
Redshift Survey: 

• Luminosity Normalization: b(L) — 0.85 + 0.15(L/L»); 

• Red Galaxy Subset Normalization: &(red) = 1.3[0.85 + 0.15(i/L»)]; and 

• Blue Galaxy Subset Normalization: fr(blue) = 0.9[0.85 + 0.15(L/Z*)]. 

Obtaining b(0) from numerical simulations of galaxy formation is also not simple. 
It is difficult to build models for galaxy populations of dark halos that can robustly 
relate the amplitude of large-scale galaxy clustering at better than the 10% level. 1^21 
This, in part, is the reason that we have a 10% limit on the accuracy of our test, 
which examines whether b = 1 at z = for a viable DE model. 

The numerical simulations of the Virgo Consortium are consistent with our 
test for DE models. Very high-res olut ion simulations have been made by the Virgo 
Consortium for the ACDM model.SH They found that for galaxies with Mb < — 17 
at z = on the largest scales, the galaxy power spectrum has the same shape as that 
of the DM, but with a slightly lower amplitude, corresponding to a bias b = 0.92. 
Samples of brighter galaxies have a bias close to b ~ 1.0. 

In this article, we related b(0) to the ACDM value of b, &a(0), using a function F. 
The bias parameter b is related to the factor F by F = \{b 2 /b\) — l\ at z = (Eq. 16). 
Thus, a maximum deviation of 10% of b 2 (z = 0) from b\(z = 0) with b\(z = 0) ~ 
1.0 implies a ~ 10% deviation of F from zero or F max = 0.1. We investigated 
DE models that make a negligible contribution to the total energy density before 
the recombination era and, thus, a negligible contribution to the density power 
spectrum normalization factor As at the recombination era. The function F has 
the important characteristic of being independent of As ■ We calculated the value of 
F numerically for several well-known DE models from the growth equation for 6p/p. 
The constraints from the Gold SNIa datak^and the condition that F — 0.1 restrict 
the values of the parameters of the linear EOS, w(a) = wq + w a (l — a), for dark 
energy. It was found that the best fit values of wq and w a are —1.86 < wq < —1.72, 
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with 1.53 < w a < 2.9. For z ~ 0.5 — 1, where w{a) is sensitive to the supernova 
data, and wo ~ — 1- 

The BWM and VMM were then studied using the factor F. We showed that these 
DE models do not satisfy our DE model test, with F = 0.1 and Q° M = 0.28 ± 0.02. 
The BWM has F = 0.26 and the VMM, F = 0.20. Thus BWM and VMM arc not 
viable DE models. 

Finally, we analyzed the SUGRA model for the above parametrized EOS with 
wo = —0.82 and w a — 0.58. F was found to be very large: -Fsugra ~ 0.38 io!o2 f° r 
Q% = 0.301q;q2, giving b 2 SUGRA (z = 0) = 1.38, which is appreciably greater than 
unity. Thus the SUGRA model is also not a viable DE model. 
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